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We discuss the use of Wilson fermions with twisted mass for simulations of QCD thermody- 
namics. As a prerequisite for a future analysis of the finite-temperature transition making use 
of automatic 6 [a) improvement, we investigate the phase structure in the space spanned by the 
hopping parameter K, the coupling j3, and the twisted mass parameter jj,. We present results for 
Nf = 2 degenerate quarks on a 16 3 x 8 lattice, for which we investigate the possibility of an Aoki 
phase existing at strong coupling and vanishing ji, as well as of a thermal phase transition at 
moderate gauge couplings and non-vanishing ji. 
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1. Introduction 



Wilson fermions with a twisted mass term evaluated at "maximal twist" offer a number of ad- 
vantages over untwisted fermions, such as the absence of exceptional configurations and automatic 
&{a) improvement. In order to see whether these features can also be successfully exploited for 
future analyses of the QCD finite temperature phase transition, we present as a prerequisite an in- 
vestigation of the phase diagram spanned by the hopping parameter K, the twisted mass parameter 
/I and the gauge coupling j8 for Nf = 2 quarks. 

The twisted mass action for two degenerate fermions consists of the standard Wilson action, 
augmented by a twisted mass term \j7iLij 5 T 3 Y- The hopping parameter K and the twist parameter 
jj. are related to the bare quark mass as 



m,. 




(1.1) 



For maximal twist we have fc = K c , in which case the quark mass is determined by jX alone. For a 
recent review of the twisted mass formulation, see [p. 

The expected phase structure is based on symmetry arguments and leading order chiral per- 
turbation theory ^ |J. It has recently been discussed in [H], a qualitative sketch is shown in 
Figure RL The phase diagram divides into two regions. The expected Aoki phase for small values 
of j6 (cf. [|]]) is in the ix = plane, since }X^0 explicitly breaks the associated parity-flavour sym- 
metry. For larger values of /3 , a conical surface is believed to mark the thermal transition between 
the hadronic and quark gluon regimes. This is based on the expectation [0] that, for fixed /3 , the 



lines of equal bare quark mass in the (fc,/i) plane given by equation ( |P| ) should correspond to 
lines of constant physics at that lattice spacing. 





3.3 3.4 3.5 3.6 3.7 3.8 3.9 



Figure 1: Speculative phase structure for twisted Figure 2: Summary of results. Note that the 

mass Wilson fermions. The Aoki phase is part of left part of the diagram shows the jX = plane 

the jX = plane and the thermal transition is ex- whereas the right part (also the enlarged insert) 

pected to be located on a conical surface; cf. [|J]. shows the thermal transition line for jX = 0.005. 

Our simulations are performed on a 16 3 x 8 lattice using an improved HMC code as described 
in [0], based on the tree-level Symanzik improved gauge action. We search for signatures for 
the Aoki phase and its possible decoupling from the deconfinement transition in the range from 
j8 = 3.40 to 3.60. For larger /3, we present results for the slice of the surface at fixed ix = 0.005 
obtained from simulations at /? € {3.75,3.775,3.8}, supplemented by data at /3 = 3.9 from a pre- 
vious investigation [Q]. Our findings, to be discussed in the following sections, are summarized in 
Figure [| 
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2. Search for the Aoki Phase 

In a first attempt to roughly estimate a region of parameter space for a potential Aoki phase, 
we performed a scan in K at several values of j8 keeping /J. = fixed. According to the Aoki 
phase scenario the masses of the charged as well as the neutral pions vanish along the border of the 
Aoki phase. Consequently, we expect the average number of conjugate-gradient iterations (Ncg) 
in the HMC process to peak there. Such a behaviour is demonstrated in the lower-left graph of 
Figure || where — log (Ncg) is shown as a function of k for different values of j8. We find that the 
minimum of — log(Ncc) shifts slightly towards lower K as j8 increases. In the upper-left panel of 
Figure ||, this is consistently reflected in the behaviour of the average plaquette, showing a rather 
rapid change when approaching the region where (Ncg) peaks. Therefore, given the data of (Ncg}> 
a possible Aoki phase should lie within the narrow interval K £ [0.16,0.20]. 

The Polyakov loop as an indicator for a thermal phase transition remains flat throughout this 
interval. A significant rise can only be found for K > 0.21, i.e. well beyond the possible Aoki 
phase region. This can be seen on the right hand side of Figure || where the real part of the 
Polyakov loop and its susceptibility are shown as functions of k. Therefore, if an Aoki phase 
were found somewhere between K = 0.16 and 0.20, this would be in contrast to the experience 
gained in previous investigations (with Wilson gauge action) ^ in which no clear separation of 
the Aoki phase from the thermal transition has been observed. As expected for the phase diagram 
at jJ. = 0, the fc value at which the Polyakov loop rises significantly shifts to lower values when j8 
is increased. 

At j3 = 3.4 we have started a more thorough investigation of the existence of a phase of broken 
parity-flavour symmetry by searching for a fc-region where the order parameter (xjfiysT^Xj/) does not 
vanish in the limit lim/^o hmy^oo ( 1//7/5 T 3 i//). Here the strength of the "external field" h is related 
to the twisted mass parameter pL by h = 2k\i. So far, we have only measured the order parameter 
as a function of K at j8 = 3.4 and the three values h = 0.01, 0.005 and 0.0025. Our current results 
for the average plaquette and the order parameter are shown in Figure ||; on the left hand side both 
observables as functions of k and on the right hand side only (l/n/sT 3 y) versus h. As above, the 
average plaquette exposes a rapid change or discontinuity as a function of K. This becomes slightly 
stronger with decreasing h. The order parameter is observed to peak around K* m 0.1825. This 
peak becomes smaller and narrower in width as h decreases. Looking at the right hand side of 
Figure ^, we get strong indications that for h — > the order parameter tends to zero for all K values 
considered. Therefore, our yet preliminary results at /3 =3.4 point to an absence of the Aoki phase 
and a possible first order phase transition, a scenario sustained by the occurence of long-lasting 
metastable states at K = 0.1830. In this context it should be mentioned that first order behaviour 
for Wilson fermions in the intermediate coupling regime has been observed before (cf. [Q| and 
references therein). We interprete our data in the same way as a remnant of the Kc(T = 0) bulk 
transition. Of course, the existence of an Aoki phase at smaller j8 is left open. In any case, our 
results for the Polyakov loop show that the deconfining transition happens at considerably larger K 
values. 
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Figure 3: Scan of the phase diagram at ji = looking for signals of an Aoki phase. Left: Average plaquette 
(top) and minus the logarithm of the average number of CG-iterations (bottom) as functions of re. Right: 
The real part of the Polyakov loop (top) and its susceptibility (bottom) also as functions of K. 
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Figure 4: Search for the Aoki -phase at j3 = 3.4 at 3 different values of h. Left: Average plaquette (top) and 
the order parameter (bottom) versus K. Right: The order parameter versus h for K < K* (top) and K > K* 
(bottom). 
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3. Thermal Transition Line 



In order to locate the thermal transition line for p. = 0.005 we fixed j8 and scanned in K, 
reaching 10000) statistics for each pairing. As observables we used the plaquette, the real 
part of the Polyakov loop and their respective susceptibilities as well as the pion norm, which 
is given by \\7l\\ = T,x^ r W( x )75¥( x )W{^)75¥{^)- The susceptibilities have been normalized as 
Xo = (^(O 2 ) — (0) 2 ^j. The critical coupling Kj for the thermal transition is then signalled by a 
peak in the susceptibilities. 

In Figure |6] we present our data for the search for the thermal transition line at p = 0.005 for 
the various values of /3. The critical Kj of the transition is decreasing with growing j3. This is 
consistent with the fact that T c (m) is an increasing function of quark mass, as found with staggered 
or untwisted Wilson fermions (for a recent review, see [|l0|]). Correspondingly, the signal in the 
Polyakov loop becomes more pronounced with increasing j3, whereas the signal in the plaquette 
becomes broader. The values we extracted for the thermal transition Kj- (/!,/?) are collected in 
Figure ||; they are compared to the values for K C (T = 0,j8), cf. [11, [TJ|, |TJ]]. As in the case of 
untwisted Wilson fermions JT4]], the distance between Kj- and K c grows with increasing j3. 
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Figure 5: Results for the thermal transition line 
Kt(P,H = 0.005) compared to the known points 
for the critical K c at T = 0. 
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As stated above, due to equation (|Ll|), for fixed lattice spacing, viz. coupling j3, there should 
be a line of equal physics in the (jc,ju) plane. This leads to the expectation of a conical thermal 
transition surface in the phase diagram as sketched in Figure |[ In order to check this expectation 
we have scanned for the transition both in K for fixed pi, as well as in p for fixed K. 

In Figure [7] we present the comparison of a scan in p that was performed at K = Kc(T = 0, j8 = 
3.9) with ^(4000) statistics with a scan in K at J3 = 3.9 and pi = 0.005 with up to ^(10000) HMC 
sweeps. The points are mapped to an effective bare quark mass using equation ( \l.l\ >. The similar 
behaviour provides preliminary evidence for the existence of at least part of the conical structure 
of the transition surface. Note that the data correspond to different /^-values. Since this varying of 
pL, i. e. of the twist angle arctan (ju /{m — m c )), introduces 0(a) effects, the data are distorted away 
from the unique curve by finite lattice spacing effects. 

In Figure |8| we sketch an approximate slice of the phase diagram in Figure [I] at a larger j8. A 
line of equal critical mass, defined by (|l~I|), for the thermal transition is sketched. The vertical line 
within the ellipse is the I st order line containing jc c (j8 , T = 0, pi = 0) as predicted by %PT for T = 
- cf. [ph. Our present and planned future investigations are illustrated as well. 
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Figure 6: Search for the thermal transition at 
ji = 0.005 and three values of j3. Left, from 
top to bottom: Average plaquette, real part of 
the Polyakov loop and pion norm. Right, from 
top to bottom: Susceptibility of plaquette and 
real part of the Polyakov loop. 
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4. Summary and Outlook 

For j6 = 3.4 and pL = we have found a very narrow range of K- values where we observe in- 
dications for a first order phase transition similar to the bulk transition seen in the zero-temperature 
case. There seems to be no Aoki phase in this range of the phase diagram. The deconfinement 
transition occurs at considerably higher fc-values. 

For larger j8 we have identified the thermal transition line at /I = 0.005, and provided some 
evidence for one side of the conical phase boundary predicted in [Q]. 

Our future work will consist of a continued investigation of the strong coupling region at /3 - 
values substantially lower than presented here as well as probing the details of the conical structure 
for larger /3 . Finally, we would like to study the thermal transition at maximal twist and to work 
towards physical quark masses and the continuum in this regime. 
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(k",ju) pairs at j3 =3.9 mapped to the bare quark the context of the expected phase structure, 
mass. 
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